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04 An overview is given of theoretical progress on self-organization at the nanoscale in reactive 

» I systems of heterogeneous catalysis observed by field emission microscopy techniques and at the 

Cu molecular scale in copolymerization processes. The results are presented in the perspective of recent 

^^ advances in nonequilibrium thermodynamics and statistical mechanics, allowing us to understand 

^^ how nanosystems driven away from equilibrium can manifest directionality and dynamical order. 
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I. INTRODUCTION 



o 

2i In macroscopic systems, self-organization arises far from equilibrium beyond critical thresholds where the macrostate 

r*^ issued from thermodynamic equilibrium becomes unstable and new macrostates emerge through bifurcations. Such 

•j bifurcations may lead to oscillatory behavior and spatial or spatio-temporal patterns, called dissipative structures 

C^ because they are maintained at the expense of free-energy sources .'^"21 Recent developments have been concerned 

c/) with the complexity of such nonequilibrium behavior, particularly, in small systems of nanometric size down to the 

^ molecular scale.E"^] The molecular structure of matter largely contributes to the complexity of natural phenomena by 

C^ the multiplicity and the variety of chemical species and their possible specific actions. Furthermore, the microscopic 

degrees of freedom manifest themselves at the nanoscale as molecular and thermal fluctuations, which requires a 

stochastic description for the thermodynamic and kinetic properties of small systems. Remarkably, great advances 

have been recently achieved, leading to a fundamental understanding of the emergence of dynamical order in fluctuating 

nonequilibrium systems, as overviewed in Section [Til 

^ These advances allow us to bridge the gap between the microscopic and macroscopic levels of description, especially, 

'— ' in reactions of heterogeneous catalysis studied by field electron and field ion microscopy techniques (FEM and FIM, 

^^ respectively) .'SHISl In such reactions, dynamical patterns are observed on metallic tips with a curvature radius of tens of 

►^ nanometers. Accordingly, the crystalline surface is multifaceted so that adsorption, desorption, reaction and transport 

rv^ processes have various speeds on different facets. Moreover, the activation barriers are significantly modified by the 

f — high electric field present under FEM or FIM conditions, as calculated by quantum electronic ah initio and density 

<3\ functional theories. '^^'i^ These non-uniform and anisotropic properties of the catalytic surface participate to the 

^^ generation of nanopatterns in the far-from-equilibrium regi mes of bistability and oscillations observed, in particular, 

(T^ during water formation on rhodium, as presented in Section [lllpSlSni 

^^ Under nonequilibrium conditions, the emergence of dynamical order is already in action at the molecular scale during 

CN copolymerization processes. Copolymers are special because they constitute the smallest physico-chemical supports of 

^yj information. Little is known about the thermodynamics and kinetics of information processing in copolymerizations 

although such reactions play an essential role in many complex systems, e.g. in biology. In this context, the recent 

advances in the thermodynamics o f stoc hastic processes are shedding light on the generation of information-rich 

copolymers, as explained in Section [ivp nEll 

jrt The purpose of this contribution is to present an overview of these recent advances on self-organization at the 

nanoscale in the perspective of future theoretical and experimental work on these topics, as discussed in the next 

sections. 



II. FUNDAMENTAL ASPECTS OF NONEQUILIBRIUM NANOSYSTEMS 

A. Structure and function of nanosystems 

There exist many different processes and systems at the nanoscale: heterogeneous catalysis on nanoparticlesjS'i^l 
electrochemical re actio ns on nanoelectrodesp^l synthetic molecular machines, ^^ ^^ single enzymes p ^ line ar and rotary 
molecular motors p2129| DNA and RNA polymerases responsible for replication and transcriptionj^IEIl or ribosomes 
performing the translation of mRNAs into proteins!^ Every nanosystem has a specific structure and acquires its 
function when driven out of equilibrium by some free-energy source. In this regard, the structure and function of a 
nanosystem can be characterized in terms of its equilibrium and nonequilibrium properties, as shown in Table IT] 

TABLE I: Comparison between the equilibrium and nonequilibrium properties of nanosystems. 



Equilibrium Nonequilibrium 

zero affinities non zero affinities 

zero mean fluxes non zero mean fluxes 

zero entropy production positive entropy production 

no free-energy supply needed free-energy supply required 

detailed balancing directionality 

3D spatial structure 4D spatio-temporal dynamics 

structure function 

A nanosystem can be in thermodynamic equilibrium with its environment at given temperature and chemical 
potentials, as it is the case for a catalytic surface in contact with a gaseous mixture at chemical equilibrium or for 
an enzyme in a solution also at chemical equilibrium. In these equilibrium systems, the ceaseless movements of 
thermal and molecular fluctuations do not need the supply of free energy to persist. In particular, the catalytic sites 
are randomly visited by adsorbates or substrates but, on average, there is no flux of matter or energy between the 
pools of reactants and products. The ratio of partial pressures or concentrations remains at its equilibrium value 
fixed by the mass action law of Guldberg and Waage. Any movement in one direction is balanced by the reversed 
movement according to the principle of detailed balancing. At equilibrium, there is no energy dissipation and no 
entropy production. From a statistical-mechanical viewpoint, the molecular architecture of the nanosystem can be 
characterized in terms of the average relative positions of its atoms but their average velocities are vanishing. In this 
respect, the nanosystem has only a 3D spatial structure at equilibrium. 

In contrast, a nanosystem in an environment containing a mixture which is not in chemical equilibrium for the 
reactions the nanosystem can catalyze will sustain non- vanishing fluxes of matter or energy. These average movements 
are driven by the free-energy sources of the environment. Energy is dissipated and entropy produced. For instance, 
a Fi-ATPase molecular motor rotates in a specific directi on if it is surrounded by a solution containing an excess of 
ATP with respect to the products of its hydrolysis .l^SESES] Therefore, nonequilibrium nanosystems acquire an average 
directionality, which can be controlled by the external nonequilibrium constraints, and they perform a 4D spatio- 
temporal dynamics, which is the expression of their function. Accordingly, the function of a nanosystem holds in 
the specific 4D dynamics that its 3D structure can developed when it is driven away from equilibrium under specific 
conditions. 

B. Out-of-equilibrium directionality of fluctuating currents 

The kinetics of nanosystems close to or far from equilibrium has recently known tremendous progress with the 
advent of time-reversal symmetry relations also called fluctuation theorems. These results flnd their ori gins in the 
study of large-deviation properties of chaotic dynamical systems sustaining transport processes of diffusion .l^^*^ Several 
versions of such relations have been obtained for systems under transient or stationary nonequilibrium conditions .I^SIHSI 
A particular version of the fluctuation theorem concerns the nonequili brium work on single molecules subjected to the 
time-dependent forces of optical tweezers or atomic force microscopy.'^Si^ For nonequilibrium systems in stationary 
states, a general fluctuation theorem has been proved for all the currents flowing across open stochastic or quantum 
systems by using microreversibility.ESHSil 

Nanosystems can be driven out of equilibrium by several independent thermodynamic forces, also called affinities.'^ 



For isothermal systems, they are defined as 



AC 

A = T^ (7=l,2,...,c) (1) 

in terms of the Gibbs free-energy differences AG^ — G^ — G^^cq supphed by the nonequihbrium environment to 
power the mean motion. They are external control parameters that depend on the concentrations or partial pressures 
of reactants and products. A nanosystem between reservoirs at different temperatures and chemical potentials is 
characterized by thermal as well as chemical affinities!^^ 

The affinities drive the cur rents flowing across the system. Examples of such currents are the reaction rates'*^ or the 
velocity of a molecular motor.E^MlAt the microscopic level of description, the instantaneous currents j(t) = {^^(i)}?,^! 
can be defined when particles cross fictitious surfaces separating reactants from products,SZI as in reaction rate theory. 
The instantaneous currents can be integrated over some time interval t to get the numbers of particles having crossed 
the fictitious surface during that time interval: AA^-y — L j~j{t') dt' . As long as the time interval t is finite, the currents 
defined as 

AA' 1 /■* 
J,^^ = -J Mt')dt' (7 = l,2,...,c) (2) 

are random variables. For given values of the different independent affinities A — {Aj}'^^i, the nanosystem reaches a 
stationary state in the long-time limit t — >■ +oo. This stationary state is described by a probability distribution Pa- 
Since the currents are fluctuating, they may take positive or negative values J = {J-y}j^i with certain probabilities 
^a(J). 

Now, we compare the probabilities of opposite values for the currents, Pa(J) and Pa(^J)- In general, these 
probabilites are different but, most remarkably, the time-reversal symmetry of the underlying microscopic dynamics 
implies that the ratio of these probabilities has a general behavior expressed by the current fluctuation theorem;^ 

^^^"^^ ~exp(A-Ji) for t ^ +(X) (3) 



^a(-J) 



This result holds for the equilibrium as well as the nonequihbrium stationary states at any value of the affinities in 
Markovian or semi-Markovian stochastic processes if the large-deviation properties of the process are well defined in 
the long-time limit .EHMl 

At equilibrium where the affinities vanish, the exponential function takes the unit value and we recover the prin- 
ciple of detailed balancing according to which the probabilities of opposite fluctuations are equal: Po(J) — Po{—-i)- 
However, out of equilibrium when the affinities do not vanish, the ratio of probabilities typically increases or decreases 
exponentially in time depending on the sign of A • J = X]-y=i ■^jJ'y- Therefore, a bias grows between the probabilities 
of opposite fluctuations and the current fluctuations soon become more probable in one particular direction. Direc- 
tionality has thus appeared in the system. This directionality is controlled by the affinities because the currents would 
flow in the opposite direction if the affinities were reversed, as expected from microreversibility. 

The current fluctuation theorem has several implications. As a consequence of Jensen's inequality, (e~^) > e^^^\ 
the thermodynamic entropy production is always non negative: 



1 diS 
fee dt 



A • (J)a > (4) 



where (J) a are the mean values of the currents in the stationary state described by the probability distribution Pa- 
Therefore, the second law of thermodynamics is the consequence of the current fluctuation theorem. Furthermore, 
this theorem allows us to generalize the Onsager reciprocity relations and the Green-Kubo formulas from the linear 
to the nonlinear response properties of the average currents with respect to the affinities. -^ -- This generalization is 
the result of the validity of the current fluctuation theorem far from equilibrium. 

In particular, these results apply to the reversible Brusselator model of oscillatory reactions.^ For fully irreversible 
reactions, in which the rates of the reversed reactions vanish, the corresponding affinities take infinite values so that the 
entropy production is also infinite, in which case the ratio of the probabilities of opposite fluctuations is consistently 
either zero or infinite. The fact that the ratio ((sl) behaves exponentially means that the reversed processes may soon 
become so rare that their probabilities are negligible and the system be in a far-from-equilibrium regime which could 
be considered as fully irreversible. 



C. Thermodynamic origins of dynamical order 

Another time-reversal symmetry relationship concerns the statistical properties of the histories or paths followed 
by a system under stroboscopic observations at some sampling time Ai. Such observations generate a sequence of 
coarse-grained states: 

W = W1W2 ■■■^n (5) 

corresponding to the successive times tj — j At with j — 1,2, ...,n. This history or path has a certain probability 
Pa('*') to happen if the system is in the stationary state corresponding to the affinities A. Because of the randomness 
of the molecular fluctuations, these path probabilities typically decrease exponentially as 

PA(a;)-PA(c^i^2---c^„)^e-"^*''^ (6) 

at a rate /ia that characterizes the temporal disor der in the process. Such a characterization concerns stochastic 
processes as well as chaotic dynamical systemsl^^Ml j^ nonequilibrium stationary states, the time-reversed path 

U>^ ^UJn-- ■ UJ2UJ1 (7) 

is expected to happen with a different probability 

PA(a;^) = PaK • • •c^2c^i) ^ e-"^*''A (8) 

decreasing at a different rate h^ now characterizing the temporal disorder of the time-reversed pathsi^The remarkable 
result is that the difference between the disorders of the time-reversed and typical paths is equal to the thermodynamic 
entropy production!^ 



1 djS 

kvt dt 



hl-hp^>Q (9) 



The second law of thermodynamics is satisfied because this difference is known in mathematics to be always non 
negative!S3 The validity of the formula (l9|) has been verified in experiments where the nonequilibrium constraints are 
imposed by fixing t he cur rents instead of the affinities, in which case the comparison should be carried out between 
Pj(a;) and P_j(a;R).lSSIS2l 

At equilibrium, detailed balancing holds so that every history and its time reversal are equiprobable, their temporal 
disorders are equal, and the entropy production vanishes. This is no longer the case away from equilibrium where 
the typical paths are more probable than their time reversals. Consequently, the time-reversal symmetry is broken at 
the statistical level of description in terms of the probability distribution Pa of the nonequilibrium stationary state. 
In this regard, the entropy production is a measure of the time asymmetry in the temporal disorders of the typical 
histories and their time reversals. As a corollary of the second law, we thus have the theorem of nonequilibrium 
temporal ordering;^" 

In nonequilibrium stationary states, the typical histories are more ordered in time than their corresponding time 
reversals in the sense that h^ < h^. 

This temporal ordering is possible out of equilibrium at the expense of the increase of phase-space disorder so that 
there is no contradiction with Boltzmann's interpretation of the second law. The result established by this theorem is 
that nonequilibrium processes can generate dynamical order, which is a key feature of biological phenomena. In partic- 
ular, the nonequilibrium ordering mechanism can generate oscillatory behavior in surface reactions or information-rich 
sequences in copolymers. 

III. HETEROGENEOUS CATALYTIC REACTIONS IN HIGH ELECTRIC FIELDS 

Dynamical order in the form of nonequilibrium patterns or oscillations can arise in heterogeneous catalysis at the 
nanoscale on metallic tips under FEM or FIM conditions.'^"!^ 

A. Surface conditions in FEM and FIM 

Since solid metals are crystalline and the radius of curvature of typical tips may reach 10-30 nm, the surface is 
multifaceted, which introduces non-uniformities as shown by Gerhard Ertl and coworkers.™ Moreover, the surface is 



subjected to high electric fields of about 10 V/nm, which have influence on the surface reactions.l^ The electrostatic 
edge effect tends to concentrate the electric field on the sharpest structures of the metallic needle. This is the case at 
the edges of the crystalline facets, which provides the atomic resolution of cryogenic FIM!^ At a larger scale, this is 
also the case at the apex of the tip where the radius of curvature of the average surface is the smallest. If the average 
shape of the tip is a paraboloid, the electric field varies as 

- ^^ (10) 



1 -I- 2— 

as a function of the radial distance r from the axis of cylindrical symmetry of the paraboloid. The quantity R denotes 
the radius of curvature at the apex where the electric field reaches its maximum value _Fo • 

This high electric field has several effects, which create the conditions of a nanoreactor localized near the apex of 
the needle under electric tension. On the one hand, the electric field polarizes the molecules in the gaseous mixture 
around the needle. Consequently, the partial pressures increase according to 

P(F)^P(0)exp^ (11) 

where a is an effective polarizability of the molecules of a given species and fee Boltzmann's constant .1^^ On the other 
hand, the electric field modifies the activation energies of the different processes taking place on the surface: 

E^{F) = E^{Q) - d^F - ^F^ + ■■■ (12) 

The dependence of the activation energies is in general a nonlinear function of the electric field F. At sufficiently 
low values of the electric field, such power expansions define the coefficients dx, a^,---, which are associated with the 
transition state in analogy to the situation in the stable states. 

Since the surface is multifaceted, its properties also depend on the orientation of every crystalline plane where 
the reactions proceed. This dependence can be expressed in terms of the corresponding Miller indices (h, fc, I) or, 
equivalently, the unit vector normal to the plane: 

n=(n.,n„n.).^=i||iL= (13) 



This anisotropy concerns, in particular, the activation energies, which can be expanded in kubic harmonics as 

for face-centered cubic crystals such as rhodium or platinum. The coefficients of this expansion can be fitted to data 
collected for different orientations. The knowledge of any activation energy for the three main crystalline planes (001), 



(Oil), and (111) determines the three first coefficients of the expansion (14 1. Finer dependences of the activation 
energy can be included with experimental data on more crystalline planes. These dependences on the various facets 
composing the tip are crucial to understand the anisotropy of the surface reactions and the nanopatterns observed 
under FEM or FIM conditions. 

These conditions are significantly different from those prevailing on flat crystalline surfaces extending over distances 
of several hundreds or thousands of micrometers. Such flat surfaces have uniform properties so that the transport 
mechanisms are diffusive and the patterns observed on flat crystalline surfaces emerge as the result of standard 
reaction-diffusion processes in uniform media. The wavelengths of such reaction-diffusion patterns are determined 
by the diffusion coefficients D and the reaction constants fcrxn- They are of the order of L ~ -^/ZJ/fcrxn ~ 100 /im, 
which is much larger than the size R ~ 20 nni of a FIM tip. Therefore, the nanopatterns observed under FEM or 
FIM conditions are not standard reaction-diffusion patterns and their understanding requires to take into account 
the non-uniform and anisotropic effects of the electric field and the underlying crystal. On non-uniform surfaces, the 
transport of adsorbates is not only driven by the gradients of coverages, but also by the surface gradients of desorption 
energy and electric field.'^^lMl 

B. Adsorption-desorption kinetics 

For diatomic molecules A2 such as dihydrogen, adsorption is dissociative. Therefore, the coverage 9a of the surface 
by the atomic species increases at the rate: 



dOA 
dt 



= 2fcaPA. (I-^a)' with fca^ , % =^ (15) 

ads •^/ZTrmAaft^BJ 



where the pressure Pa2 is given by Eq. (11), S° is the initial sticking coefficient at zero coverage, and Og is some 
reference area. Accordingly, desorption is dissociative and proceeds at the thermally activated rate 

d0A 



dt 



= -2fcd^i with fcd^fc^e-'^^^'"'-^^^*) (16) 

dcs 



and /3 — (fceT)^^. The desorption energy E^ depends on the electric field (lOl, the local crystalline orientation (13l, 
as well as the coverage itself if lateral interactions play a role in desorption. 

Experimental data are available for the adsorption and desorption of hydrogen on rhodium 

H2 (gas) + 20 (ad) ^^ 2 H (ad) (17) 

Its desorption energy takes values below one electron- Volt so that the mean hydrogen coverage 6'h is low above 400 K 
where bistability and oscillations are observed.'i^"221 

C. Surface oxides of rhodium 

The behavior is more complex for oxygen on rhodium. On the one hand, the dissociative adsorption of oxy gen 
involves a precursor state, which explains the dependence of the sticking coefficient on the oxygen coverage.'^^'^ 

On the other hand, oxygen forms surface oxide trilayers on rhodium: 0(ad)-Rh-0(sub). Their properties have 
been systematically investigated, in particular, with DFT calculations .ESHzSiBy their stoichiometry Rh02, the surface 
oxides differ from the bulk oxide Rh203. Oxygen vacancies can thus exist in either the outer or the inner oxygen 
layer. Accordingly, the structure of a partially formed surface oxide can be characterized in terms of the occupancies 
by oxygen atoms of the adsorption and subsurface sites, Oq and ^s respectively. 

These different features of the interaction of oxygen with rhodium are described by the following kinetic schemefl^"2ni 

adsorption and desorption: O2 (gas) + surface _^ O2 (pre) + surface 

k o 

dissociation and recombination: O2 (pre) +20 (ad) ^20 (ad) 

fcdO 

k 

oxidation and reduction of Rh: O (ad) + (sub) ^ (ad) + O (sub) 

fcrod 

(18) 

The surface oxide trilayer tends to inhibit oxygen adsorption, which is taken into account by the precursor 
const ant ;^Sj'20i 

with parameters fitted to data on the sticking coefhcient) ^^ * ^^ * Besides, the desorotion rate constants on the three main 
surface orientations are fitted to temperature-programmed desorption spectra.^ISHSSl 

For rhodium in equilibrium with gaseous dioxygen, the adsorbate and subsurface occupancies satisfy 



^O fcrcd ^s / ^aO ^aO „ /„„x 



1 — Oq fcox 1 - ^s V ^dO fcdO 

A phase transition between a metallic surface with adsorbed oxygen and a surface oxide trilayer with possible vacancies 
occurs at about 6*3 — 0.5. This condition allows us to determine the ratio of the oxidation and reduction rates of a 
rhodium layer using the results of DFT calculations on the three main surface orientations.'iSHini 

Moreover, the dependence of the activation energies on the electric field has been evaluated using DFT calculations,'^^ 
which shows that a positive electric field tends to promote the oxidation of rhodium. 

D. The H2-02/Rh system 

Water formation is catalyzed on rhodium in contact with a gaseous mixture of dihydrogen and dioxygen. If water 
vapor is not supplied in the mixture, its partial pressure is vanishing. Under these conditions, the corresponding 
affinity M is infinite and the overall reaction proceeds in a fully irreversible regime. 



For the reaction of water formation on rhodium 



2 H (ad) + O (ad) ^30 (ad) + H2O (gas) 



(21) 



the rate constant is taken as 

fcr = fc°e"'^(^'"'*'-^+^"''H+'4?''o) (22) 

with coefficients A^ and A^ expressing the change of the reaction rate with the hydrogen and oxygen coverages, 
respectivelyP^lf22l 



Kinetic equations 



Combining together the different processes involving hydrogen and oxygen, the kinetic equations are given by^SHSQI 



dt 
dt 

dt 



2 fcaH Ph2 
2 



- 2 fcdH 6'h - 2 k, 9u 9o - div Jr 



1 + Ke^^ 

koxOo (1 - 
with the coverage of empty sites 00 

Jh = -Dt 

with the effective energy potential 



k,.oKP, 



O2 C0 



fcdO ^O 1 ^ ^' 



■Cox CO 



(1 - 9^) + fcrod ( 



fcr On ec 



- fc. 



edi 



(23) 
(24) 

(25) 



1 



9o and K = ^ao/^do- The current density of hydrogen takes the form 



(1 - 9o)Ven + 0hV0o + dMl -Oh- ^o) 



knT 



^dH(r)-ddHi^(r) + -aH,F(r)' 



est 



(26) 



(27) 



The mobility of oxygen is negligible. In contrast, the mobility of hydrogen is very high with a typical diffusion time 
tdiS = R^ /D-H ~ 10~^ s, which is many orders of magnitude shorter than the time scales of the other kinetic processes 
at 500 K. Therefore, the coverage of atomic hydrogen quickly reaches the quasi-equilibrium distribution 



?H(r,t) = 



1 



)(r,0 



„/J[(7H(r)-MH(i)] 



1 



(28) 



such that the hydrogen current density is vanishing: Jh — 0. This coverage varies in space because the effective 
potential (27) depends on the electric field (Il0| and on the normal unit vector (13), which is uniquely determined 



by the position r on the surface. Therefore, the hydrogen coverage already forms nanopatterns, which reflect the 
anisotropy of the underlying crystal and the electric field variation. The chemical potential /in is uniform on the 
surface but slowly varies in time because the population of hydrogen atoms on the surface changes due to adsorption, 
desorption and the reaction of hydrogen with oxygen. Using multiscale analysis, it has been possible to show that 
this nonequilibrium chemical potential evolves in time according to 



dt 



= knT- 



Efacets^[(i-^o)atfe + fea, 



'too 



Efacets^fe(l-^H-eC 



(29) 



where w = {1 — 9o) ^, while the time derivatives dt9ii and dfOo are given by the kinetic Eqs. (23) and (24) 



^F 



This kinetic model reproduces very well the nonlinear dynamics of the system and the spatial dependence of the 
observed nanopatterns. 



2. Bistability 

Bistability manifests itself under variation of hydrogen pressure for fixed oxygen pressure, as shown in Fig. [T] as 
a function of temperature.1SH221 At low (resp. high) hydrogen pressure, the surface is covered with oxygen (resp. 
hydrogen). A domain of hysteresis appears in between where the two states coexist. The bistability domain depends 
on the applied electric field. The higher the field, the broader the coexisting region in the bifurcation diagram of 
Fig. IT] As temperature increases, so does the water formation rate, leading to the reduction of the oxide and the 
coverage of the tip by hydrogen. Ultimately, the bistability domain disappears above 550 K. 




400 450 500 

Temperature (K) 



FIG. 1: The bistability diagram showing hysteresis for P02 = 5 x 10~* Pa during catalytic water formation on rhodium at (a) 
11.0 V/nm and (b) 12.3 V/nm. The circles and stars indicate the experimental pressures for which the structural transformation 
occurred when decreasing and increasing the hydrogen pressure, respectively. The area in betw een marks the coexistence region 
of bistability. The full lines are the corresponding results for the kinetic model ruled by Eqs. (23l-(29l.'^^'^ 



3. Oscillations 



The kinetic model also explains the oscillatory behavior observed in this system in FIM experiments (see Figs. [2] 
and pi) .'^SJEni The period of oscillations is about 40 s. In the model, this period is mainly determined by the rate con- 
stants of rhodium oxidation and reduction when oxygen reacts with the first rhodium layer. The feedback mechanism 
at the origin of the oscillations involves, in particular, the formation of surface oxide and its inhibition of further 
oxygen adsorption, as taken into account with the rate constant (19). 




FIG. 2: Series of FIM micrographs covering the complete oscillatory cycle as well as the corres ponding time evolution of the 
subsurface oxygen distribution on a logarithmic scale as obtained within the kinetic modeL'^^'^ The temperature, electric field 
and partial pressures of oxygen are T = 550 K, Fq = 12 V/nm, P02 = 2 x 10~^ Pa, respectively. On the other hand, the 
hydrogen pressure is Ph2 = 2 x 10""^ Pa in the FIM experiments and 4 x 10~^ Pa in the simulation of the kinetic model. For 
the subsurface site occupation, the white areas indicate a high site occupation value while the dark areas indicate a low site 
occupation value. 



Starting from a quasi metallic surface in Figs. [2k & [2H, an oxide layer invades the topmost plane and grows along 
the {011} facets forming a nanometric cross-like structure seen in Figs. l2]3 & [2]3. The oxide front spreads to cover 
finally the whole visible surface area in Figs. [2j: & [2f . This is associated with a decrease of the overall brightness 
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FIG. 3; (a) The time evolution of the oxygen coverage (solid line) and the oxygen subsurface occupation (dashed line) at the 
(001) plane in the kinetic model in the oscillatory regime, (b) Corresponding oscillations of the hydrogen coverage at the (001) 
plane, (c) Experimental total brightness during the oscillations.^^ The conditions are the same as in Fig. M The arrows indicate 
the transition from a metallic rhodium field emitter tip to one that is invaded by subsurface oxygen. 

in Fig. |3j:. The oscillation cycle is closed by a sudden reduction of the surface oxide from the outskirts towards the 
top, with a considerable increase of the brightness. During the cycle shown in Fig. [31 rhodium is alternatively covered 
by the surface oxide when 0s — 1 and 9o — 0.5 and, thereafter, by hydrogen at a low coverage since the temperature 
is 550 K. 



E. Self-organization at the nanoscale 



In summary, the nonlinear dynamics and the patterns observed in field emission microscopy are determined by the 
tip geometry, the anisotropy from the underlying crystal, and the electric field. The anisotropy can be described by 
giving to the energy barriers their dependence on the crystalline orientation for the many facets composing the tip 
of the field emission microscope, in particular, with the systematic expansions ( 14 ) into kubic harmonics for cubic 



crystals. Another important feature is the ultrafast mobility of hydrogen, which remains in the quasi-equilibrium 
distribution (28) slowly evolving because of the other kinetic processes. The kinetic model is built on the basis of 



experimental data about adsorption and desorption of hydrogen and oxygen as revealed by temperature-programmed 
desorption spectra and the recent studies on the Rh02 rhodium surface oxides.^^"^ By taking into account all these 
different aspects, the kinetic model provides a comprehensive understanding of the bistability, the oscillations, and 
the nanopatterns observed in FIM experiments.'i^"2Sl 

Chemical nanoclocks have been o bserved under field emission microscopy conditions in several reactions besides 
water formation on rhodiumj ^^ l ^^^^^^ Surprisingly, rhythmic behavior is possible at the nanoscale of 10-30 nm in the 
population dynamics of the different species reacting on the surface. Indeed, this area may contain up to about ten 
thousand adsorbates, which is already much larger than the minimum size of a few hundred molecules required to 
sustain correlated oscillations .^^Zl This behavior is an example of dynamical order, which can manifest itself out of 
equilibrium as a corollary of the second law. 



IV. COPOLYMERIZATION PROCESSES 



A. Information processing at the molecular scale 



If the history of a nonequilibrium system can be recorded on a spatial support of information, the theorem of 
nonequilibrium temporal ordering**^ suggests that dynamical order may generate regular information sequences, which 
is not possible at equilibrium. 



10 

At the molecular scale, natural supports of information are given by random copolymers where information is coded 
in the covalent bonds. This is the idea of Schrodinger's aperiodic crystal.'^ Random copolymers exist in chemical 
and biological systems. Examples are styrene-butadiene rubber, proteins, RNA, and DNA, this latter playing the role 
of information support in biology. Accordingly, dynamical aspects of information are involved in copolymcrization 
processes where fundamental connections with nonequilibrium thermodynamics have been recently discovered ."^l' ■^-2, 

B. Thermodynamics of free copolymerization 

The stochastic growth of a single copolymer proceeds by attachment and detachment of monomers {m} continuously 
supplied by the surrounding solution, which is assumed to be sufficiently large to play the role of a reservoir where 
the concentrations of the monomers are kept constant: 

to = mim2 • ■ -rni ^ w = mim2 ■ ■ ■ mimi+i (30) 

The probability Pt{^) to find the monomer sequence lo of length I = \u)\ at the time t is ruled by the master equation 

dPt{uj) 



dt 



J2 [■Pt(w') W{uj'\uj) - Pt{uj) W{uj\uj')] (31) 



where the coefficients M^(a;|a;') denote the rates of the transitions cu — ^ ui' . If attachment and detachment processes 
are slower than the equilibration of the chain with its environment, the transition rates satisfy the conditions of local 
detailed balancing 

in terms of the Gibbs free energy G{uj) of a single copolymer chain uj in the solution at the temperature T. The 
enthalpy Hioj) and the entropy S{lo) of the copolymer chain lo can similarly be defined and they are related together 
by G{uj) = H{u}) — TS{uj). At a given time t, the system may be in different sequences and different configurations 
so that the total entropy has two contributions: 

St=Y, PtHS{uj) - fcB ^ Ptiio) InPt(w) (33) 

The first one is due to the statistical average of the phase-space disorder S{lo) of the individual copolymer chains uj 
and the second is due to the probability distribution itself among the different possible sequences uj existing at the 

current time t. 

In the regime of steady growth,E2180| this probability is supposed to be factorized as 

Ptiuj) ^ fiiiu;) X ptil) (34) 

into a stationary statistical distribution fii{i-u) describing the arbitrarily long sequence which is left behind, multiplied 
by the time-dependent probability pt{l) of the length I selected in the sequence. In this regime, the mean growth 
velocity is constant and given by w = d{l)t/dt where {l)t = '^iPt{l) x I is the mean length of the copolymer at time 
t. Mean values per monomer can be defined as 



for entropy s, enthalpy h, and Gibbs free energy g — h — Ts 



\im\y ^Jil{Lu)X{oJ) (35) 



In these circumstances, the total entropy (33 1 can be calculated and shown to vary in time as 

dSt _ dcS_ d]S_ , , 

dt ~ dt dt ^ ' 

due to the entropy exchange between the copolymer and its surrounding: 

dcS h 
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FIG. 4: Space-time plot of the stochastic growth of a copolymer composed of two monomers, as simulated by Gillespie's 
algorithm with the parameter values: fe+i = 1, fc_i = 10~^, [1] = 10~^, fc+2 = 2, fc_2 — 2 x 10^^, and [2] = 5 x 10~*, in 
the growth regime by entropic effect!^ Under these conditions, the fraction of monomers 1 is p = 0.618, the growth velocity 
V = 6.17 X 10""* (i.e., 0.18 monomer per reactive event), the free-energy driving force e — -~g/{k-B,T) = —0.265, the Shannon 
disorder D — 0.665, the affinity A = £-\- D — 0.400, and the entropy production -^ — Av — 2.47 x 10~* in units where /sb ~ 1. 



and the entropy production: 






kjiAvX) 



(38) 



which is always non negative according to the second law of thermodynamics. The entropy production is expressed 
in terms of the afRnit}^^!] 



A = - 



knT 



-D (polymer) 



(39) 



which involves, on the one hand, the free energy per monomer g and, on the other hand, the Shannon disorder per 
monomer in the sequence composing the copolymer: 



Z)(polymer) = lim — - > /^;(w) ln/i;(w) > 



(40) 



A prediction of this result is that a copolymer can grow by the entropic effect of disorder Z? > in an adverse free- 
energy landscape as long as the affinity (391 is positivci^^^This is illustrated in Figs.Hand^for free copolymerization 
with two monomers. The concentration [2] of the second monomer is kept fixed while the concentration [1] of the 
first one is varied. The growth velocity, as well as the affinity (39) and the entropy production (38 1, all vanish at 
equilibrium for [1] = [l]oq, which does not coincide with the concentration [1] = [l]o where the free-energy driving 
force £ = —g/{k-QT) is vanishing. Therefore, the growth is possible for intermediate values of the concentration 
[l]eq < [1] < [l]o in the entropic growth regime, preceding the growth regime driven by free energy when e > for 
[1] > [l]o. At equilibrium, the Shannon disorder (40) reaches its maximum value _D = ln2 and decreases away from 
equilibrium, as seen in Fig. [5^. 



C. Thermodynamics of copolymerization with a template 



A similar result holds for copolymerizations with a template. If the sequence of the template a is characterized 
by the statistical distribution vi{a), the Shannon conditional disorder of the copy uj with respect to the template is 
defined aJ^H 



D(copy|template) = lim — 7 / z/i(a) /i;(cLi|a) ln/i;(aj|a) > 



(41) 
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FIG. 5: Comparison between simulation (dots) and theory (lines) for the growth of a copolymer cornposed of two monomers 
with the parameter values: fc+i — 1, k-i = 10~^, fc+2 = 2, fc_2 = 2 x 10~^, and [2] = 5 x 10~''M' Several characteristic 
quantities are depicted versus the concentration [1] of monomers 1: (a) the fraction p of monomers 1 in the copolymer (circles) 
and the Shannon disorder D (squares); (b) the growth velocity v (triangles), the free-energy driving force e (open circles), the 
affinity A = s + D (open squares), and the entropy production -^ = Av (crossed squares). The entropy production vanishes 
at the equilibrium concentration [l]cq = 5 x 10"* together with the velocity and the affinity. However, the free-energy driving 
force vanishes at the larger concentration [l]o = 1.30453 x 10"'^. The regime of growth by entropy effect exists between these 
two values of the concentration. The velocity and the affinity are positive for [1] > [l]eq = 5 x 10~^ and negative below (not 
shown). The free-energy driving force is positive for [1] > [l]o ~ 1.30453 x 10""^ and negative below (not shown). 



and the mutual information between the copy and the template asSZl 

/(copy, template) = _D(copy) — Z? (copy | template) > (42) 

where the Shannon disorder of t he c opy is defined as in Eq. ( |40[ ). In this framework, the thermodynamic entropy 
production is again given by Eq. ( 38 ) but with the affinity^ 



A — — h I?(copy|template) 

fcei 



knT 



D{copy) — /(copy, template) 



(43) 



which establishes quantitatively a fundamental link between information and thermodynamics at the molecular scale. 



D. The case of DNA replication 



The previous results apply to the different living copolymerization processes and, in particular, to DNA replication. 
In this case, the subunits of the polymers are the four nucleotides N = A, T, C, or G and the monomers the 
corresponding nucleoside triphosphates NTP. Assuming that no free energy difference exists between correct and 
incorrect chains, the copolymerization process can be simulated by Gillespie's algorithm as a function of the driving 
force e = —g/{kBT)^ The results are depicted in Fig. p^ which shows the percentage of replication errors as well as 
the mutual information between the copy and the template. The error percentage is maximum at equilibrium and it 
decreases as the growth is pushed away from equilibrium. Similarly, the mutual information vanishes at equilibrium 
and saturates at /max = 1-337 nats for high enough values of the driving force. As in the case of free copolymerization, 
a transition occurs between the regime of growth by entropic effect for eoq < £ < and the growth driven by free 
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FIG. 6: Stochastic simulation^ of DNA replication by polymerase Pol 7 of human mitochondrial DNA from GenBank^'^ using 
known data on the kinetic constants of Watson-Crick pairing.^ The reversed kinetic constants are taken as k-mn = k+mn e"*^. 
(a) Percentage of errors in DNA replication versus the driving force e. (b) Mutual information between the copied DNA strand 
and the original one versus the driving force. The arrow points to the equilibrium value of the driving force: Seq — — ln4 = 
-1.38629. 

energy for e > 0. At equilibrium, information transmission is not possible between the template and the copy. Fidelity 
in the copying process becomes possible if enough free energy is supplied during copolymerization. 

The existence of growth by entropic effect could be experimentally investigated in chemical or biological copolymer- 
izations. In polymer science, methods have not yet been much developed to perform the synthesis and sequencing of 
copolymers for the information they may support. However, such methods are already well developed for DNA and 
under development for single-molecule DNA or RNA sequencing.^^ ■^•^ These methods could be used to test experi- 
mentally the predictions of copolymerization thermodynamics by varying NTP and pyrophosphate concentrations to 
approach the regime near equilibrium where the mutation rate increases. 



CONCLUSIONS AND PERSPECTIVES 



Nowadays, self-organization has been studied for different phenomena from the macroscopic world down to the 
molecular scale. At the macroscale, self-organization emerges far from equilibrium beyond bifurcations leading to 
bistability or oscillatory behavior.'i"^! However, on smaller and smaller scales, the time evolution of physico-chemical 
systems is more and more affected by thermal and molecular fluctuations which are the manifestation of the microscopic 
degrees of freedom. In the framework of the theory of stochastic processes, the probability to find the system in some 
coarse-grained state is ruled by a master equation. The macroscopic description in terms of deterministic kinetic 
equations is only obtained in the large-system limit, in which bifurcations emerge between nonequilibrium macrostates. 
This emergence concerns in particular oscillations which only become correlated if the system is large enough;^^ 

This is the case for the nanoclocks of heterogeneous catalysis observed by field emission microscopy techniques 
such as the reaction of water formation on rhodium in high electric field under FIM conditions .I^Hini Jn the present 
contribution, methods are described for the modeling of these nonequilibrium processes at the nanoscale on highly 
non-uniform and anisotropic surfaces. The underlying crystalline structure determines the reactivity of the surface 
and contributes to the formation of the observed nanopatterns, which cannot be interpreted as standard reaction- 
diffusion patterns. In spite of their nanometric size, such systems are large enough to undergo self-organization and 
manifest bistability and rhythmic behavior. 

This dynamical order is the result of directionality which is induced away from equilibrium by the external con- 
straints and free-energy sources characterized by the empowering affinities. If the principle of detailed balancing holds 
at equilibrium, this is no longer the case out of equilibrium so that the fluctuating currents are biased and acquire a 
directionality, which is expressed by the current fluctuation theorem! "^^ 1 This theorem has for consequence the non- 
negativity of the entropy production in accordance with the second law of thermodynamics, as well as generaliza tions 
of the Onsager reciprocity relations and Green-Kubo formula from the linear to the nonlinear response properties.HSEIl 
Furthermore, the entropy production appears as a measure of time asymmetry in the temporal disorder of the typical 
histories of a system. Out of equilibriu m, th e typical histories are more probable than their time reversals, which 
appears as a corollary of the second law.ESinSl These new results transcend the known formulation of nonequilibrium 
thermodynamics, explicitly showing how the second law finds its origin in the breaking of time-reversal symmetry at 
the mesoscopic level of description in the theory of nonequilibrium systems. 

These results show that the second law of thermodynamics already governs self-organization at the molecular 
scale. Thermodynamics can be applied to the stochastic growth of a single copolymer, which is the support of 
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information encoded in its sequence. Away from equilibrium, dynamical order enables information processing during 
copolymerizations, which is not possible at equilibrium. The statement by Manfred Eigen that "information cannot 
originate in a system that is at equilibrium'^SSl is rigorously proved in this framework. The thermodynamics of 
copolymerization shows that, indeed, fundamental connections exist between information and thermodynamics.'^i''221 
In particular, the growth of a copolymer can be driven either by the free energy of the attachment of new monomers 
or by the entropic effect of disorder in the grown sequence. 

These considerations open new perspectives to understand the dynamical aspects of information in biology. During 
copolymerization processes with a template as it is the case for replication, transcription or translation in biological 
systems, information is transmitted although errors may occur due to molecular fluctuations, which are sources of 
mutations. Metabolism and self-reproduction, which are the two main features of biological systems, turn out to 
be linked in a fundamental way since information processing is constrained by energy dissipation during copolymer- 
izations. Moreover, the error threshold for the emergence of quasi-species in the hypercycle theory by Eigen and 
Schustej^ could be induced at the molecular scale by the transition towards high fidelity replication beyond the 
threshold at zero free energy per monomer between the two growth regimes.'^ In this way, prebiotic chemistry could 
be more closely linked to the first steps of biological evolution. 

Acknowledgements 

This research has been financially supported by the F.R.S.-FNRS, the "Communautc frangaise de Belgique" (con- 
tract "Actions de Recherche Concertees" No. 04/09-312), and the Belgian Federal Government (lAP project "NOSY"). 



4 



I. Prigogine, Introduction to Thermodynamics of Irreversible Processes (Wiley, New York, 1967). 
^ G. Nicolis and I. Prigogine, Self- Organization in Nonequilihrium Systems (Wiley, New York, 1977). 
3 R. Imbihl and G. Ertl, Chem. Rev. 95, 697 (1995). 

G. Nicolis and I. Prigogine, Exploring Complexity (W. H. Freeman and Company, New York, 1989). 

G. Nicolis and C. Nicolis, Foundations of Complex Systems (World Scientific, New Jersey, 2007). 

R. Feistel and W. Ebeling, Physics of Self-Organization and Evolution (Wiley- VCH, Weinheim, 2011). 

P. Gaspard, in: G. Radons, B. Rumpf, and H. G. Schuster, Editors, Nonlinear Dynamics of Nanosystems (Wiley- VCH, 

Weinheim, 2010) pp. 1-74. 

M. F. H. van Tol, A. Gielbert, and B. E. Nieuwenhuys, Catal. Lett. 16, 297 (1992). 

V. Gorodetskii, J. Lauterbach, H.-H. Rotermund, J. H. Block, and G. Ertl, Nature 370, 276 (1994). 

N. Ernst, G. Bozdech, V. Gorodetskii, H.-J. Kreuzer, R. L. C. Wang, and J. H. Block, Surf. Set. 318, L1211 (1994). 

C. Voss and N. Kruse, Appl. Surf Sci. 94/95, 186 (1996). 

N. Kruse and T. Visart de Bocarme, in: G. Ertl, H. Knozinger, F. Schiith, and J. Weitkamp, Editors, Handbook of Hetero- 
geneous Catalysis, 2nd edition (Wiley- VCH, Weinheim, 2008) pp. 870-895. 

T. Visart de Bocarme, N. Kruse, P. Gaspard, and H. J. Kreuzer, J. Chem. Phys. 125, 054704 (2006). 

J.-S. McEwen and P. Gaspard, J. Chem. Phys. 125, 214707 (2006). 

J.-S. McEwen, P. Gaspard, F. Mittendorfer, T. Visart de Bocarme, and N. Kruse, Chem. Phys. Lett. 452, 133 (2008). 

T. Visart de Bocarme, T. Bar, and N. Kruse, Ultramicroscopy 89, 75 (2001). 

T. Visart de Bocarme, G. Beketov, and N. Kruse, Surf Interface Anal. 36, 522 (2004). 

J.-S. McEwen, P. Gaspard, T. Visart de Bocarme, and N. Kruse, Proc. Natl. Acad. Sci. USA 106, 3006 (2009). 
^^ J.-S. McEwen, P. Gaspard, T. Visart de Bocarme, and N. Kruse, J. Phys. Chem. C 113, 17045 (2009). 
^° J.-S. McEwen, P. Gaspard, T. Visart de Bocarme, and N. Kruse, Surf Sci. 604, 1353 (2010). 
^^ D. Andrieux and P. Gaspard, Proc. Natl. Acad. Sci. USA 105, 9516 (2008). 
^2 C. Jarzynski, Proc. Natl. Acad. Sci. USA 105, 9451 (2008). 
2^ D. Andrieux and P. Gaspard, J. Chem. Phys. 130, 014901 (2009). 
2" V. Garcia-Morales and K. Krischer, Proc. Natl. Acad. Set. USA 107, 4528 (2010). 
^'' J. -P. Sauvage and P. Gaspard, Editors, From Non-Covalent Assemblies to Molecular Machines (Wiley- VCH, Weinheim, 

2011). 
^® P. Lussis, T. Svaldo-Lanero, A. Bertocco, C.-A. Fustin, D. A. Leigh, and A.-S. Duwez, Nature Nanotech. 6, 553 (2011). 
^"^ W. Min, L. Jiang, J. Yu, S. C. Kou, H. Qian, and X. S. Xie, Nana Lett. 5, 2373 (2005). 
^* T. Funatsu, Y. Harada, M. Tokunaga, K. Saito, and T. Yanagida, Nature 374, 555 (1995). 
2^ H. Noji, R. Yasuda, M. Yoshida, and K. Kinosita Jr., Nature 386, 299 (1997). 
^° K. Adelman, A. La Porta, T. J. Santangelo, J. T. Lis, J. W. Roberts, and M. D. Wang, Proc. Natl Acad. Sci. USA 99, 

13538 (2002). 
^^ W. J. Greenleaf and S. M. Block, Science 313, 801 (2006). 
^2 J. Eid et al. Science 323, 133 (2009). 



18 



15 

^^ S. Uemura, C. Echeverria Aitken, J. Korlach, B. A. Flusberg, S. W. Turner, and J. D. Puglisi, Nature 464, 1012 (2010). 

^'' D. Andrieux and P. Gaspard, Phys. Rev. E 74, 011906 (2006). 

^^ E. Gerritsma and P. Gaspard, Biophys. Rev. Lett. 5, 163 (2010). 

^•^ P. Gaspard and G. Nicolis, Phys. Rev. Lett. 65, 1693 (1990). 

^■^ P. Gaspard, New J. Phys. 7, 77 (2005). 

^* D. J. Evans, E. G. D. Cohen, and G. P. Morriss, Phys. Rev. Lett. 71, 2401 (1993). 

^8 G. Gallavotti and E. G. D. Cohen, Phys. Rev. Lett. 74, 2694 (1995). 

■'o G. E. Crooks, Phys. Rev. S 60, 2721 (1999). 

■'^ J. L. Lebowitz and H. Spohn, J. Stat. Phys. 95, 333 (1999). 

^^ M. Esposito, U. Harbola, and S. Mukamel, Rev. Mod. Phys. 81, 1665 (2009). 

■^^ C. Jarzynski, Annu. Rev. Condens. Matter Phys. 2, 329 (2011). 

'*'' D. Cohin, F. Ritort, C. Jarzynski, S. B. Smith, I. Tinoco Jr., and C. Bustamante, Nature 437, 231 (2005). 

'^^ D. Andrieux and P. Gaspard, J. Chem. Phys. 121, 6167 (2004). 

"^ D. Andrieux and P. Gaspard, J. Stat. Mech.: Th. Exp. POlOll (2006). 

''■^ D. Andrieux and P. Gaspard, J. Stat. Phys. 127, 107 (2007). 

''* D. Andrieux and P. Gaspard, J. Stat. Mech.: Th. Exp. P02006 (2007). 

"^ D. Andrieux and P. Gaspard, J. Stat. Mech.: Th. Exp. P11007 (2008). 

^° D. Andrieux and P. Gaspard, J. Stat. Mech.: Th. Exp. P02057 (2009). 

■'^ D. Andrieux, P. Gaspard, T. Monnai, and S. Tasaki, New J. Phys. 11, 043014 (2009). 

^^ T. De Bonder and P. Van Rysselberghe, Affinity (Stanford University Press, Menlo Park CA, 1936). 

^^ D. Andrieux and P. Gaspard, J. Chem. Phys. 128, 154506 (2008). 

'"^ P. Gaspard and X.-J. Wang, Phys. Rep. 235, 291 (1993). 

^^ P. Gaspard, Chaos, Scattering and Statistical Mechanics (Cambridge University Press, Cambridge UK, 1998). 

"^^ P. Gaspard, J. Stat. Phys. 117, 599 (2004). 

^^ T. M. Cover and J. A. Thomas, Elements of Information Theory, 2nd edition (Wiley, Hoboken, 2006). 

^^ D. Andrieux, P. Gaspard, S. CiHberto, N. Gamier, S. Joubaud, and A. Petrosyan, Phys. Rev. Lett. 98, 150601 (2007). 

''^ D. Andrieux, P. Gaspard, S. CiHberto, N. Garnier, S. Joubaud, and A. Petrosyan, J. Stat. Mech.: Th. Exp. P01002 (2008). 

^° P. Gaspard, C R. Physique 8, 598 (2007). 

^^ E. W. Miiller and T. T. Tsong, Field Ion Microscopy: Principles and Applications (Elsevier, New York, 1969). 

®^ A. S. Mikhailov and G. Ertl, Chem. Phys. Lett. 238, 104 (1995). 

®^ M. Hildebrand, A. S. Mikhailov, and G. Ertl, Phys. Rev. £58, 5483 (1998). 

^^ Y. De Decker, A. Marbach, M. Hinz, S. Giinther, M. Kiskinova, A. S. Mikhailov, and R. Imbihl, Phys. Rev. Lett. 92, 198305 

(2004). 
^'^ A. Garcia Cantii Ros, J.-S. McEwen, and P. Gaspard, Phys. Rev. £'83, 021604 (2011). 
®® E. Schwarz, J. Lenz, H. Wohlgemuth, and K. Christmann, Vacuum 41, 167 (1990). 
^"^ A. N. Salanov and V. I. Savchenko, Surf. Sci. 296, 393 (1993). 

^^ M. V. Ganduglia-Pirovano, K. Reuter, and M. Scheffler, Phys. Rev. B 65, 245426 (2002). 
®^ J. Gustafson, A. Mikkelsen, M. Borg, E. Lundgren, L. Kohler, G. Kresse, M. Schmid, P. Varga, J. Yuhara, X. Torrelles, C. 

Quiros, and J. N. Andersen, Phys. Rev. Lett. 92, 126102 (2004). 
™ C. Africh, F. Esch, W. X. Li, M. Corso, B. Hammer, R. Rosei, and G. Comelh, Phys. Rev. Lett. 93, 126104 (2004). 
^^ L. Kohler, G. Kresse, M. Schmid, E. Lundgren, J. Gustafson, A. Mikkelsen, M. Borg, J. Yuhara, J. N. Andersen, M. 

Marsman, and P. Varga, Phys. Rev. Lett. 93, 266103 (2004). 
'^^ J. Gustafson, A. Mikkelsen, M. Borg, J. N. Andersen, E. Lundgren, C. Klein, W. Hofer, M. Schmid, P. Varga, L. Kohler, 

G. Kresse, N. Kasper, A. Stierle, and H. Dosch, Phys. Rev. B 71, 115442 (2005). 
''^ C. Dri, C. Africh, F. Esch, G. ComeUi, O. Dubay, L. Kohler, F. Mittendorfer, G. Kresse, P. Dubin, and M. Kiskinova, J. 

Chem. Phys. 125, 094701 (2006). 
^^ E. Lundgren, A. Mikkelsen, J. N. Andersen, G. Kresse, M. Schmid, and P. Varga, J. Phys.: Condens. Matter 18, R481 

(2006). 
■^^ F. Mittendorfer, J. Phys.: Condens. Matter 22, 393001 (2010). 

''^ J.-S. McEwen, P. Gaspard, Y. De Decker, C. Barroo, T. Visart de Bocarme, and N. Kruse, Langmmr 26, 16381 (2010). 
■^■^ P. Gaspard, J. Chem. Phys. 117, 8905 (2002). 

^* E. Schrodinger, What is Life? (Cambridge University Press, Cambridge UK, 1944). 
™ B. D. Coleman and T. G. Fox, J. Polym. Set. A 1, 3183 (1963). 
*° B. D. Coleman and T. G. Fox, J. Chem. Phys. 38, 1065 (1963). 

*^ Homo sapiens DNA mitochondrion, 16569 bp, locus AC 000021, version CI: 1 15315570, |htt p : //www . ncbi . nlm . nih . gov 
**^ H. Lee and K. Johnson, J. Biol. Chem. 281, 36236 (2006). 

*^ M. Eigen, Steps towards Life: A Perspective on Evolution (Oxford University Press, Oxford, 1992). 
** M. Eigen and P. Schuster, Naturwissenschaften 64, 541 (1977); ibid. 65, 7 (1978); ibid. 65, 341 (1978). 
*^ H.-J. Woo and A. Wallqvist, Phys. Rev. Lett. 106, 060601 (2011). 



